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SOLUTION OF THE INVERSE SPECTRAL PROBLEM 
FOR A CONVOLUTION INTEGRO-DIFFERENTIAL OPERATOR 
WITH ROBIN BOUNDARY CONDITIONS 


S.A. Buterin0 and A.E. Cheque Rivero@ 

Abstract. The operator of double differentiation on a finite interval with Robin boundary 
conditions perturbed by the composition of a Volterra convolution operator and the differentia¬ 
tion one is considered. We study the inverse problem of recovering the convolution kernel along 
with a coefficient of the boundary conditions from the spectrum. We prove the uniqueness 
theorem and that the standard asymptotics is a necessary and sufficient condition for an 
arbitrary sequence of complex numbers to be the spectrum of such an operator. A constructive 
procedure for solving the inverse problem is given. 
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1. INTRODUCTION 

Consider the boundary value problem L = L{M, h, H) of the form 

iy-.=—y" + f M {x — t)y'(t) dt = Xy, 0 < x < tt, (1) 

Jo 

U{y) := y\0) - hy{0) = 0, V{y) := y'{7i) + Hy^n) = 0, (2) 

where A is the spectral parameter, M{x) is a complex-valued function, {tt—x)M{x) GL2(0,vr) 
and h, H ^ C. 

We study an inverse spectral problem for L. Inverse problems of spectral analysis consist 
in recovering operators from given their spectral characteristics. Such problems often appear 
in mathematics, mechanics, physics, electronics, geophysics, meteorology and other branches of 
natural sciences and engineering. The greatest success in the inverse problem theory has been 
achieved for the Sturm-Liouville operator iiy := —y" + q{x)y (see, e.g., [1-4]) and afterwards 
for higher-order differential operators [5-7]. For example, it is known that the potential q{x) 
can be uniquely determined by specifying the spectra of two boundary value problems for 
equation iiy = Xy with one common boundary condition. 

For integro-differential and other classes of nonlocal operators inverse problems are more 
difficult for investigation, and the classical methods either are not applicable to them or 
require essential modifications (see [4, 8-19] and the references therein). In [9] a perturbation of 
the Sturm-Liouville operator with Dirichlet boundary conditions by the Volterra convolution 
operator was considered. It was proven that the specification of only the spectrum uniquely 
determines the convolution component. Moreover, developing the idea of Borg’s method a 
constructive procedure for solving this inverse problem was obtained along with the local 
solvability and stability. In [15] the global solvability was proved by reducing this inverse 
problem to solving the so-called main nonlinear integral equation, which was solved globally. 
Earlier by a particular case of this approach the analogous results were obtained for the 
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operator (1) with Dirichlet boundary conditions [12], Here we study the case of Robin boundary 
conditions. A short version of this preprint is to appear in [20]. 

Let (p{x, A) be a solution of equation ( 1 ) satisfying the initial conditions 

(^(0,A) = 1, ¥p'(0,A) = h. (3) 

Clearly, the eigenvalues of the problem L with account of multiplicity coincide with the zeros 
of the function 

A(A):=R(^(a;,A)), (4) 

which is called the characteristic function of L. By the wellknown method (see, e.g., [4]) 
involving Rouche’s theorem one can prove that the spectrum of L consists of inhnitely many 
eigenvalues A^, /c > 0. Moreover, the following theorem holds. 

Theorem 1. The spectrum {Afc}fc>o has the form 

= {k + Kk)'^, {Kk} G l2- (5) 


As compared with the Dirichlet boundary conditions the Robin ones (2) bring additional 
difficulties in studying the inverse problem for L. First, we consider the following problem. 

Inverse Problem 1. Given {Afc}fc>o and h, H, hnd M{x). 

For this inverse problem we prove the following uniqueness theorem. 

Theorem 2. The specification of the spectrum {Afc}fc>o uniquely determines the function 
M{x), provided that the coefficients h, H are known a priori. 

We note that Inverse Problem 1 is overdetermined. Indeed, the spectrum possesses also 
some information on the coefficients of the boundary conditions. In particular, we prove that if 
h = 0, then along with M{x) the coefficient H is also determined, i.e. a uniqueness theorem 
holds for the following problem. 

Inverse Problem 2. Given {Afc}fc>o and h = 0, find H and M{x). 

Moreover, the following theorem holds. 

Theorem 3. For arbitrary complex numbers Xk, k > 0, of the form (5) there exists a 
unique (up to values on a set of measure zero) function M{x), (tt — x)M{x) G L 2 ( 0 , 7 r), and 
a unique number H E C, such that {Afc}fc>o is the spectrum of the problem L{M,0, H). 

Thus, the asymptotics (5) is a necessary and sufficient condition for the solvability of Inverse 
Problem 2. The importance of the assumption h = 0 is explained in Remark 1 (see Section 4). 
We leave open whether the analogous criterium can be obtained assuming only that h is known 
a priori but h 7 ^ 0, or symmetrically: H is given while h is unknown. 

In the next section we derive the main nonlinear integral equation of the inverse problem 
and prove the global solvability of this nonlinear equation. In Section 3 we prove some auxiliary 
assertions along with Theorem 1. In Section 4 we give the proof of Theorems 2 and 3, which is 
constructive, and provide algorithms for solving the inverse problems (Algorithms 1 and 2 ). 

2. MAIN NONLINEAR INTEGRAL EQUATION 

Let the functions C{x, A), S{x, A) be solutions of equation ( 1 ) satisfying the initial conditions 
G(0, A) = A'(0, A) = 1, C"(0, A) = ^(0, A) = 0. 
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Thus, according to (3), (4) we have A) = C{x, A) + hS{x, A) and 

A(A) = C'{7i, A) + hS'{7r, A) + HC{7i, A) + hHS{7r, A). 

Let = X. The following representation is wellknown; 

S(,, A) = + r P(^, tfXlfilxJlM, 0 < X < X, 

P Jo P 

where P{x,t) is the kernel of the transformation operator. In [12] it was shown that 


U=1 




where 


N*^{x) = N{x), N<'^+^\x) = N*N*'^{x)= / N{x - t)N*^ {t) dt, z/> 1, 


( 6 ) 


(7) 


( 8 ) 


and the function iV(a;), (tt — a;)iV(a;) G ^ 2 ( 0 , tt), is connected with M{x) by the relation 

M{x) = 2N{x) — [ dt f N{t — r)iV(r) dr, 0 < x < tt. 

Jo Jo 

The following lemma gives further representations, which we use in the sequel. 
Lemma 1. The following representations hold: 

rx 

S{x, X) = K{x, x) — p / K{x,t) sin p{x — t) dt, 

Jo 

pX 

S'{x, A) = R{x, x) — p R{x, t) sin p{x — f) dt, 

Jo 

/ X 

Q{x, t) sin p{x — t) dt. 


where 


C'{x, A) = —p sin px — p P{x, t) sin p{x — t) dt. 


K{x,t)=t + / {t — t)P{x,t) dr, 

Jo 

R{x,t) = l+ f P{x,T)dT+ [ (t — t)Px{x,t) dr, 
Jo Jo 

Q{x,t) = l+ [ P{x — t + t,t) dr. 

Jo 


(9) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


Proof. The integration by parts in (7) gives 


S{x,X)= f (^ + f P{x,t) dr^ cos p{x — t) dt. 


( 17 ) 
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Integrating (17) by parts we arrive at (10) and (14). Then differentiating (10) with respect to 
X we get 

d r r 

S'{x, X) = —K{x,x) — p / KJx,t) sin p(x — t) dt — X / K{x,t) cos p{x — t) dt. 
dx Jq Jq 

Integrating by parts in the last term we arrive at 

d /■* 

S'(x, X) = —K(x,x) — p / (Kx(x,t) + Kt(x,t)) sinp(x — t) dt. 
dx Jq 

Taking (14) and (15) into account we get (11). 

Further, we have 

Cix, A) = 1 - A [ S{t, A) dt. 

Jo 

Indeed, put 

/ X px 

M{x — t)S{t,X) dt, M(a;) := 1 — A / S{t,X)dt. 


(18) 


Since n(0) = m( 0) and v\x) = u'{x), we have v{x) = u{x). On the other hand, since 

u\x) = —XS{x,X), u"{x) = —XS'{x,X), 


we get 

iu = x(^S'{x, ~ J ~ i)S{t, A) dt^ = Xv{x) = Xu{x). 

Taking into account m( 0) = 0(0, A) and m'( 0) = O'(0, A) we arrive at (18). 
Differentiating (18) and substituting (7) therein we get (13). 

Finally, substituting (7) into (18) we get 


0(a;,A) 


1 


p / sin ptdt — p 




t) sin pr dr. 


Changing the order of integration yields 


0(a;,A) 




X + t) dr ] sin p{x 


t) dt, 


which gives (12) and (16). 

The next lemma is a direct corollary of formulae (6), (10)-(13). 
Lemma 2. The characteristic function has the form 


□ 


A(A) = —psinpTT + q; + p / w{x) sinpxdx, w{x) E L 2 { 0 ,n). 


Here 

a = hR{TT, Tr,N) + H + hHK{TT, tt; N), 

—w{'K — x) = P{tt, x; N) + hR{'K, x; N) + HQiji, x; N) + hHKiji, x; N), 
where we add the argument ” N ” in order to indicate the dependence on N{x). 


(19) 

( 20 ) 
( 21 ) 



5 


The relation (21) can be considered as a nonlinear equation with respect to N{x). We 
call it main nonlinear integral equation of the inverse problem. The main equation (21) can be 
rewritten in the explicit form. Indeed, by virtue of (8), (14)-(16) we have 

CJO / \i/ 

P(7r,x;M)^y'y—^N“'(x), (22) 

U=1 


1 px "I px 

i?(7r, x; M) = 1 + — / [n — tyN*^ {t) dt+ — t){7r — ty~^N*'^{t) dt, (23) 

u=i do do 

CXD / \i/ pX 

QU, li M) = 1 + 5^ VZj£L / N-U) dt, (24) 

u=l do 

oo ^ 

K (tt, x; M) = X + — 

V = 1 

Substituting (22)-(25) into (21) we get 

OO 

/(^) = {yyx)N*''{x) + 

u=l 

where f{x) = —w{ti — x) — h — H — hHx, ^uix) = (tt — xY/nl and 

'^Yxy) = - xY + h(7r - tY~^ ^TT - t + (x - t)(z/ + Tr(7r - t))j j. 

Theorem 4. For eaeh funetion f{x) G L2(0,7r) and any complex numbers h, H equation 
(26) has a unique solution N{x), (tt — x)N{x) G L2(0,7r). 

Proof. By virtue of Theorem 4 in [10], equation (26) has a unique solution N{x), which 
belongs to L2(0,T) for each T G (0,7r). Following [10] we represent this solution in the form 

N{x) = NYx) + NYx), 

where Ni{x) E L2{0,7i) and iV 2 (a;) = 0 on (0,7r/2). Then we have 

N*Yx) = NY{x) + * N2{x), u>2. 

Substituting this into (26) we arrive at 


^Yx,t)N*Yt) dt), 


(26) 


{x — t){7i — tYN*''{t) dt. 


(25) 


TT 


f{x)-jai{x) = {7r-x)N2{x)+ / A{x,t) {tt - t)N 2 {t) dt, —<x<ti, 

J-k /2 d 


where 


AH( 2 ^) = ^ U)u{x)NY{x) + / '^u{x,t)NY{t) dt). 


U=1 

00 


= - — xi'^ijxY) + y^^iy(YYx)N(^'' ^\x-t)+ / '^Yx,t)N(^'' ^^(r-t)(irjj 


u=2 


are square-integrable functions. Hence (tt — x)N 2 {x) G ^2(0,71). 


□ 
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3. AUXILIARY ASSERTIONS 


In this section for convenience of the reader we prove Theorem 1 and further auxiliary 
assertions, which we use in Section 4 for solving the inverse problems. We note that Theorem 1 
is a direct corollary of the following assertion. 

Lemma 3. Any function of the form (19) has infinitely many zeros Xk, k > 0, having the 
asymptotics (5). 

Proof. For p E Gs '■= {p '■ \p — k\ > 6, k E Z}, 5 > 0, we have | sinp7r| > exp(|Imp|7r), 
Cs > 0. Therefore, for sufficiently large |p|, p E Gs, the following inequality holds 


|psinp7r| > 


i*7r 

a + p w{x) sin pxdx 

Jo 


According to Rouche’s theorem there are exactly A + 1 zeros Xk, k = 0,N, of the function 
A(A) lying inside the contour rAr = {A: |A| = (A+1/2)^} for sufficiently large N. Moreover, 
for each <5 > 0 there exists ks such that for \k\ > ks there is exactly one zero pk of the 
function A(p^) inside the contour 7fc((5) = {p ■ \p — k\ = 5}. Thus, pk = y/Xf = k + Hk, 
where Kk = o(l). Substituting this into (19) we get {nk} E h and (5) is proved. □ 

Analogously to Theorem 1.1.4 in [4] using Hadamard’s factorization theorem one can obtain 
the following assertion. 

Lemma 4. The function A(A) is uniquely determined by its zeros by the formula 


°° \ _ \ 

A(A) = ^(A„-A)n^ 


(27) 


Proof. It follows from (19) that A(A) is entire in A of order 1/2, and consequently by 
Hadamard‘s factorization theorem, A(A) is uniquely determined up to a multiplicative constant 
by its zeros: 



where s > 0 is the multiplicity of the eigenvalue A = 0. Consider the function 


Ao(A) := -psinpTT = -Avr]^ G ~ Ja 


k=l 


X 


Then 




^ TT 

TT AA A; 2 -A Afc , ^ 

A^.=0 


Ao(A) 


Xk-k^ 

k^-X 


where Uk 


On the other hand, taking (5) and (19) into account we calculate 


lim ^^ = 1, 

A—>■—oo 


Xk — kf 


A—A—oo^-J- \ — 


k^-X 


= 1 , 


e, kEn, 

1, A; = 0. 



and hence 
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Substituting this into (28) we arrive at (27). □ 

Thus, according to Lemmas 3 and 4 any function of the form (19) is uniquely determined by 
formula (27) from its zeros, which, in turn, have the asymptotics (5). By the standard approach 
(see, e.g.. Lemma 3.3 in [12]) one can prove the following inverse assertion. 

Lemma 5. (i) Let arbitrary complex numbers Xk, k>0, of the form (5) be given. Then the 
function A(A) determined by (27) has the form (19) with certain a E C and complex-valued 
function w{x) G L2(0,7r). 

Proof. Once again using the function 



Ao(A) := —psmpn = — 

ttA 

OO 

n 

k=l 

k^-\ 

’ 

(29) 

from (27) we get 








A(A) = Ao(A)F(A), 

F(A) = 

A 

-A 

A 

0 TT Xk — X 

J-J- k^-X’ 

k=l 

(30) 

Put 








(3k := A(P2), 

9 ■ 

t>k ■— — 

k 

a 

k>l. 

(31) 

By virtue of (29), 

(30), we have 






(Jk = 1 

(-l)‘-"‘^(Ao-A.=)(A,- 

■ k‘^)dk, 

dk 

= 

°° X 

n Xp. 

(32) 





U 

= 1, 


Hence, according 

to (5), (31), (32) we get 







9k = “ k‘^){2nkk + K,l)dk - 

Since the sequence {dk} is bounded, we have {6k} G h- Determine the function w{x) G 
7/2(0, tt) such that 


and denote 


Consider the function 


9k = w{x)smkxdx, k > 1, 

Jo 

/ TT 

w{x) sin px dx. 


(33) 


5l(p) _ _ sinpvr g(n) q, 

S{p) := + F(p2) + 


- 1 , 


sm pTT 


sm pTT 


p sm pTT 


which is, by (31) and (33) after removing the singularities, entire in p. Let us show that 
|F(p^)| < for p E Gs = {p '. \p — k\ > 5, k E Z}, <5 > 0. Denote n^k '■= —K-k, k E N. 

Then 

9 \ 

Kk 


F(/)=n (i+ 


p- 


|fc|GN 


k — p 


For a hxed 5 > 0 choose N such that < S/2 for k > N. We have for p E Gs 

.2 


F(p^) = exp(77jv(p)) ^ n (^ + 17^)’ 




0 <|A:|<Af 


(34) 



where 


Hn{p) = E In (l + ^ 

|A:|>Ar 




k — n 


Kk 


k — p ^ z/ + 1 \k — p 

v=Q 

For p E Gs using the Cauchy-Bunyakovsky inequality we arrive at 

1/2 


\Hn{p)\<c( y : n- 


|fc|>V 


\k-p\^^ 


and consequently H]sr{p) is bounded in Gs- Thus, from (34) it follows that \F{p‘^)\ < Gs for 
p G Gs- Moreover, it is easy to see that F{p^) — )■ 1 as \Imp\ — )■ cx). Using the maximum 
modulus principle we conclude that S{p) is bounded and consequently, according to Liouville’s 
theorem, it is constant. Since S{p) —)■ 0 for |/mp| — )■ cxd, S'(p) = 0 and formula (19) holds. 
□ 


4. SOLUTION OF THE INVERSE PROBLEM 

By virtue of (19), (27), the number a can be found by the formula 


a = 


ttAo JJ 


k=l 




(35) 


According to (19) the function w{x) can be reconstructed as the Fourier series 


2 ^{k) 71 — X 

w[x) = — 2_^ —:— sm kx — a- 


TT ' k 
k=l 


TT 


(36) 


Proof of Theorem 2. According to Lemmas 2 and 4 the specihcation of the spectrum uniquely 
determines the function w{x). Since h and H are given, by virtue of Theorem 4 the function 
N[x) is a unique solution of the main equation (26). Hence, the function M{x) is determined 
uniquely and can be reconstructed by formula (9). □ 

Let the spectrum {Afc}A:>o of a problem L(M, h, H) along with the numbers h, H be 
given. Then the function M{x) can be constructed by the following algorithm. 

Algorithm 1. (i) Construct the function w{x) by formulae (27), (35) and (36); 

(ii) hnd N{x) by solving the main equation (26); 

(hi) construct M{x) by formula (9). 

Proof of Theorem 3. Using the given numbers \k we construct the function A(A) by 
formula (27). According to Lemma 5 it has the representation (19) with certain number a 
and function w{x) G 7/2(0, tt). Let N{x) be the solution of the main equation (26) with this 
function w{x) and h = 0, P[ = a. Determine the function M{x), {ti — x)M{x) G 7/2(0,tt), 
by formula (9) and consider the corresponding boundary value problem L = L{M, 0, 77). It is 
easy to see that the constructed function A(A) is the characteristic function of this problem L. 
Thus, the spectrum of the latter coincides with {Afc}fc>o. According to (9), (19) and (27) the 
uniqueness of M{x) follows from the uniqueness of the solution of equation (26). □ 

Let a sequence {Afc}fc>o of the form (5) be given. According to Theorem 3 there exists 
a unique boundary value problem L{M, 0,77) with the spectrum {Afc}fc>o, which can be 
constructed by the following algorithm. 
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Algorithm 2. (i) Having calculated the number a by formula (35), put 

h := 0, H := a; 

and construct the function w{x) by formulae (27), (36); 

(ii) find N{x) by solving the main equation (26); 

(hi) construct M{x) by formula (9). 

Remark 1. The importance of the assumption h = 0 can be seen from formula (20). 
Indeed, since N{x) is unknown, the number a being, in turn, uniquely determined by the 
spectrum determines H only if h = 0. 
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